NONLINEAR DISPERSION EQUATIONS: 
SMOOTH DEFORMATIONS, COMPACTONS, 
AND EXTENSIONS TO HIGHER ORDERS 



VICTOR A. GALAKTIONOV 

Abstract. The third-order nonlinear dispersion PDE, as the key model, 
(0.1) u t = (uu x ) xx in R x R + , 

is studied. Two Riemann's problems for (jO.ip with initial data Sz^(x) — =Fsignx, create 
the shock (u(x,t) = S-(x)) and smooth rarefaction (for data 5+) waves, [18] . 

The concept of "^-entropy" solutions (a"5-entropy test") and others are developed 
for distinguishing shock and rarefaction waves by using stable smooth <5-deformations of 
discontinuous solutions. These are analogous to entropy solutions for scalar conservation 
laws such as Ut + uu x = 0, developed by Oleinik and Kruzhkov (in M. N ) in the 1950-60s. 
The Rosenau-Hyman K(2, 2) (compacton) equation 

u t = {uu x ) xx + Auu Xl 

which has a special importance for applications, is studied. Compactons as compactly 
supported travelling wave solutions are shown to pass the <5-entropy test. Shock and 
rarefaction waves are discussed for other NDEs such as 

u t = (u 2 u x ) xx , u tt = (uu x ) xx , u tt = uu x , u tt t = (uu x ) xx , u t = (uu x ) xxxxxx , etc. 



Dedicated to the memory of Professors O.A. Oleinik and S.N. Kruzhkov 

1. Introduction: nonlinear dispersion PDEs and main results 

1.1. NDEs: nonlinear dispersion equations in application and general PDE 
theory. The present paper continues the study began in [IB] of odd-order nonlinear 
dispersion (or dispersive) PDEs (NDEs). The canonical model is the third-order quadratic 
NDE (the NDE-3) 

(1.1) u t = A(u) = (uu x ) xx = uu xxx + 3u x u xx in Ex(0,T), T > 0. 
Posing for ( 11.11) the Cauchy problem includes locally integrable initial data 

(1.2) u(x, 0) = u (x) in R. 

Frequently, we assume that u is bounded and compactly supported. We will also deal with 
the initial-boundary values problem in (-L, L) x IR + with Dirichlet boundary conditions. 
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Main applications concerning NDEs can be found in [T4"l [18] ; see also [T9j Ch. 4], so 
that we do not discuss these issues in detail. However, we need to stress the attention 
of the Reader to the compacton phenomena, which were not properly treated in the 
mathematical literature. 

Compact patterns and NDEs. These are known for the Rosenau-Hyman (RH) equation 

(1.3) u t = (u 2 ) xxx + (u 2 ) x , 

which is the K(2, 2) equation from the general K(m, n) family of the NDEs: 

(1.4) u t = {u n ) xxx + {u m ) x (u>0). 

References on physical applications of such NDEs are available in [181 § 1] and in [T9J 
§ 4.2]. We will check entropy properties of compactons for various NDEs of this type. 

Further applied compacton-like models are discussed in [18J. A standard definition of 
weak solutions for (11.11) is also presented there, so that we are in a position to explain our 
main targets concerning entropy-like theory of shocks. 

1.2. Plan of the paper: entropy theory (a test) via smooth deformations and 
compactons. As in [TB], we begin with discussion of some auxiliary properties of the 
NDE-3. 

Smoothing for the NDE-3. Firstly, we recall that the smoothing phenomena and results 
for sufficiently regular solutions of linear and nonlinear third-order PDEs are well know 
from the 1980-90s. For instance, infinite C°°-smoothing results were proved in [2] for a 
general linear equation of the form 

(1.5) u t + a(x,t)u xxx = (a(x,t)>c>0), 
and in [3J for the corresponding fully nonlinear PDE 

(1-6) u t + f(u xxx ,u xx ,u x ,u,x,t) = ( f Uxxx > c > 0) . 

Namely, for a class of such equations, it is shown that, for data with minimal regularity 
and sufficient decay at infinity, there exists a unique solution u(x,t) e C°° for arbitrarily 
small t > 0. Similar smoothing local in time results for unique solutions are available for 

(1.7) ut + f{D 3 u,D 2 u,Du,u,x,y,t) = in R 2 x E + ; 

see [25] and further references therein. Concerning unique continuation and continuous 
dependence properties, see [6] and references therein, and [36] for various estimates. 

The NDE: a conservation law in H^ 1 . Writing (11.11) as (see details in [THl § 1.4]) 

(1.8) (-Dl^ut + uu x = 0, 

yields the first a priori uniform bound for data Uq G if _1 (R). Namely, multiplying (11.81) 
by u in L 2 gives the conservation law 

(1.9) 5sK*)IIIr-i = KOII*-i = <*>=IMIff-i forallt>0. 

Main results . In the present paper, we propose some concepts for developing adequate 
mathematics of NDEs with shocks, which will be concluded in Section [5] by revealing 
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connections with other classes of nonlinear degenerate PDEs. It turns out that some 
NDE concepts has definite reliable common roots and can be put into the framework of 
much better developed theory of quasilinear parabolic equations. We restrict our attention 
to a key demand, how to distinguish the shock and rarefaction waves, and this is done by 
developing the so-called "^-entropy test" on solutions via smooth deformations. General 
uniqueness-entropy theory for NDEs such as (11.11) and others is shown to be illusive |15j . 
Concerning the simple canonical model (11. II) . we do the following: 

• (i) Reviewing local existence and uniqueness theory for the NDE (II. ip and, on its 
basis, developing an 5-entropy test for distinguishing shock and rarefaction waves. 

For the RH equation such as (II. 3p . we prove that: 

• (ii) Rosenau's compacton solutions are both 5-entropy and G-admissible. 

Some of related questions and results were previously discussed in a more applied and 
formal fashion in [T3J § 7] and [THl Ch. 4]. 

1.3. On extensions and other nonlinear dispersion models. The developed con- 
cepts cover a wide range of various NDEs. First of all, we should mention that the fact 
that (II .ip is degenerate at u = and hence admits compactly supported solutions (which 
is an interesting pleasant feature) makes the analysis of 5-entropy solutions and shocks 
much harder. However, shock waves exist for other non-degenerated NDEs with analytic 
coefficients. For instance, we study entropy shocks for the NDE with infinite propagation, 

(1.10) ut = ((i + ti a K) M . 

All our further NDEs admit analogous non-degenerate versions admitting shock and rar- 
efaction waves, but no finite propagation and interfaces in the Cauchy problem. 
Another related to (II .ip model to be discussed is the cubic fully divergent NDE 

(1.11) u t = {u 2 u x ) xx = | (u 3 ) xxx (the conservation law analogy is u t + u 2 u x = 0). 

We study (11.111) instead of less physically motivated "quadratic" model u t = (\u\u x ) xx 
that exhibits similar properties of shocks and rarefaction waves. 

The results on 5-entropy solutions and similarity patterns can be extended (Section E]) 
to truly quadratic non-fully divergent NDEs such as 

(1.12) u t = (uu xx ) x = uu xxx + u x u xx , 

which we call the NDE-(2,1), where 2 and 1 stand for the number of the internal and 
external derivatives in this differential form. Notice that a standard concept of weak 
solutions hardly applies to (11.121) . so that the shock S_(x) is not a weak solution. In 
order to underline once more the fact that being weak is not a necessary demand, we 
consider a formal fully nonlinear NDE 

(1.13) \utput = {uu x ) xx , where 7 > -1. 

For 7 = 0, this gives the original equation (II. ip . Obviously, for 7 7^ 0, (11.131) does not 
admit any weak formulation. Nevertheless, we show that (11.131) admits blow-up formation 
of shocks of SL-type. 
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In Section [6], we discuss the shock formation mechanism for higher-order in time NDEs, 

(1.14) u tt = (uu x ) xx and u ttt = (uu x ) xx . 

Several principal features remain the same for higher-order NDEs such as the quadratic 
fifth-order NDE (NDE-5) 

(1.15) u t = -{uu x ) xxxx or, in general, u t = (-l) m+1 D 2 x m {uu x ), m > 1; 

see Section [3 These are conservation laws in H~ 2 , or H~ m . The mathematics of particular 
similarity solutions with shocks is developed in similar lines but technically becomes more 
involved, so we have to catch the similarity profiles numerically. 

We also claim that some concepts such as smooth (^-deformation and others, developed 
for models in ID can be adapted to the iV-dimensional NDEs. In particular, the basic 
NDE (TUU) in R N takes the form 



(1.16) M< = A( M J^) = f ^Au 2 in R N x{0,T). 

2. Conservation laws: smooth 5- deformations define entropy solutions 

From now on, being sufficiently informed about formation of crucial shock and other 
singularities in the NDEs, we will start to investigate the general questions on existence 
and uniqueness of entropy weak solutions of (II. ip . As usual, we begin our discussion by 
stressing attention to key analogies with classic theory of first-order conservation laws 
such as Euler's equation from gas dynamics 

(2.1) u t + uu x = in Ix R + . 

Entropy theory for such first-order PDEs was created by Oleinik [3U [32] and Kruzhkov 
[23] (x G M. N ) in the 1950-60s; see details on the history, main results, and modern 
developments in the well-known monographs [U [51 [34]. Thus, we now apply smooth 
5-deformation concepts to these simpler PDEs considered now in Qi = R x (0, 1). 

2.1. Preliminaries: entropy inequalities and solutions for conservation laws. It 

is known from the 1950's that the Cauchy problem for general scalar conservation laws 
admits a unique entropy solution. We refer to first complete results by Oleinik (obtained 
in 1954-56), who introduced entropy conditions in ID and proved existence and uniqueness 
results (see survey [31]), and by Kruzhkov (1970) [23], who developed general non-local 
theory of entropy solutions in R^. In the general case, one of Oleinik's local entropy 
condition has the form [31] p. 106] 

(2 _ 2) u {x ^-u(, 2 ,t) < K{xuX2jt) forall x^GR, te[0,l], 

where K is a continuous function for t G [0,1]. Oleinik's local condition E (Entropy) 
introduced in [32], for the model equation (12.11) corresponds to the well-known principle 
of non- increasing entropy from gas dynamics, 

(2.3) u{x + ,t) < u{x~,t) in Qi = Rx (0,1], 

with strict inequality on lines of discontinuity, [311 p. 101]. 



Kruzhkov's entropy condition [23J on solutions u G L°°(Qi) of (12. 1 p takes the form of 
the non-local inequality 

(2.4) \u-k\ t + \ [sign(u - k)(u 2 - k 2 )] x <0 in V'{Q X ) for any feel. 

This inequality is understood in the sense of distributions meaning that the sign < is 
preserved after multiplying the inequality by any smooth compactly supported cut-off 
function ip 6 C^(Qi), <p > 0, and integrating by parts. See clear presentation of these 
ideas in Taylor [23 p. 401]. Oleinik's and Kruzhkov's approaches are known to coincide 
in the ID geometry. Both entropy conditions generate a semigroup of contractions in L 1 , 
so that if u and v are two solutions of fl2.ll) . then 

(2.5) £K*)-«(*)IIl*<0. 

It is key that the unique entropy solution is constructed by the parabolic ^-approximation 

(2.6) u e : u t + uu x = eu xx (e > 0). 

Multiplying (12. 6p by any smooth monotone increasing function E{u) (an approximation 
of sign (u — k) for any 1:6 1) yields on integration by parts the correct sign: 

(2.7) JJ eu xx E{u) =-efJ E'(u)(u x ) 2 < 0. 

Hence, as e — > 0, this gives the necessary sign as in (12 .4p . 

The obvious advantage of the conservation law (12.11) is that, for smooth initial data 
(II .2p . the unique local continuous solution is obtained by method of characteristics and 
is given by the corresponding algebraic equation 

(2.8) dt = ^ =^ u(x,t) =u (x-u(x,t)t) for all £e[0,At), 

where At < 1 is the first moment of time when a shock of the type S^(x) (this type is 
guaranteed by (12.31) ) occurs at some point or many points. 

Thus, for t > At, it is necessary to apply the entropy inequalities to select good 
(entropy) solutions. Using this, and bearing in mind that entropy solutions are continuous 
relative initial data (in L , say), we propose the following construction which is fully based 
on algebraic relations (12. 8p : 

2.2. Conservation laws: 5-stable = entropy solutions. It is the obvious well-known 
and, nevertheless, crucial observation that, by the characteristic mechanism (12.81) . 

(2.9) non-entropy shocks of the shape S + cannot appear evolutionary. 
Indeed, differentiating (12.81) in x yields 

(2 10) u x{X,L) — 1+u > Q ( x _ u ( x> t)ty, t>u Wid>i 

u' > =^ no blow-up of u x ( "gradient catastrophe" ) occurs. 

Recalling the necessary evolution property in (I2.10p , given a small 5 > and a bounded 
(say, for simplicity, in L 1 and in L°°) solution u(x, t) of the Cauchy problem (12. ip . (II. 2p . we 
construct its 5 -deformation given explicitly by the characteristic method (12.81) as follows: 
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(i) we perform a smooth 5-deformation of initial data Uq G L 1 H L°° by introducing a 
suitable C l function u s(x) such that 

(2.11) f\u Q -u 0S \<5. 

By uis(x,t) we denote the unique local solution of the Cauchy problem with data Uqs, 
so that by (12.81) . continuous function Uis(x, t) is defined algebraically on the maximal 
interval t G [to, ti(5)), where we denote to = and ti(5) = A^. It is important that, here 
and later on, smooth deformations are performed in a small neighbourhood of possible 
discontinuities only leaving the rest of smooth profiles untouchable, so that these evolve 
along the characteristics, as usual. 

Actually, this emphasizes the obvious fact that the shocks (on a set of zero measure) 
occur as a result of nonlinear interaction of the areas with continuous solutions, which 
hence cannot be connected without discontinuities. 

(ii) Since at t = A 1( 5 a shock of type 5"_ (or possibly infinitely many shocks) is supposed 
to occur, since otherwise we continue the algebraic procedure, we perform another suitable 
5-deformation of the "data" uis(x, Ais) to get a unique continuous solution U2s(x,t) on 
the maximal interval t G [ti(5), t 2 (#)), with t 2 (<5) = Au + A 2< 5, etc. 

(k) With suitable choices of each 5-deformations of "data" at the moments t = tj(S), 
when Ujs(x, t) has a shock for j = 1, 2, there exists a t k (5) > 1 for some finite k = k(5), 
where k(5) — > +oo as 5 — ► 0. It is easy to see that, for bounded solutions, k(5) is 
always finite. A contradictions is obtained while assuming that tj(5) ^i<lasj—>-oo 
for arbitrarily small 5 > meaning a kind of "complete blow-up" that is impossible for 
conservation laws obeying the Maximum Principle. 

This gives us a global 5 -deformation in R x [0,1] of the solution u(x, t), which is a 
discontinuous orbit denoted by 

(2.12) u s (x,t) = {u jS (x,t) for t G tj(8)), j = 1, 2, .., k(6)}. 

Recall that the whole orbit (12 . 1 2[) has been constructed by the algebraic characteristic 
calculus using (12.81) only. Finally, by an arbitrary smooth 5 -deformation, we will mean 
the function (12.121) constructed by any sufficiently refined finite partition {tj(6)} of [0, 1], 
without reaching a shock of SL-type at some or all intermediate points t = tj(5). 

We next say that, given a solution u(x,t), it is stable relative smooth deformations, or 
simply 5-stable (S eformation- stable) , if for any e > 0, there exists 5 = 5(e) > such that, 
for any finite 5-deformation of u given by (12.121) . 

(2.13) ff\u-u 5 \<e. 
Then we have the following simple observation: 

Proposition 2.1. Let under given hypothesis, a weak solution u(x,t) of the Cauchy prob- 
lem (\2.1\) , (1 1 . jgp be 5 -stable. Then it is entropy. 

Indeed, if u(x,t) is not entropy, then there exists G (0, 1] such that u(x, t*) does not 
satisfy (12. 3p . i.e., this profile has a finite non-entropy shock of the type S+ at some point 
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G KL Since those shocks cannot be reproduced with arbitrary accuracy e in L l by the 
characteristic system (I2.8p . any 5-deformation u s at t = £* must stay Eq > away from 
u(x,t*) for arbitrarily small 5 > 0. 

Of course, this construction does not play a role for conservation laws with well- 
developed entropy theory, which establishes existence of a semigroup of /^-contractions of 
entropy solutions. Obviously, this strong contractivity property guarantees also unique- 
ness of ^-entropy solutions. The situation is different for the NDEs: 

3. On 5-entropy solutions (a test) of the NDE 

Thus, we are going to develop and discuss some aspects of entropy solutions for (11. ip . 
without using the idea of vanishing, e — > 0, viscosity as in [13, § 7] 

(3.1) u t = (uu x ) xx - eu xxxx in Rxl + (e > 0). 

A direct verification that the ^-approximation (13.11) yields as e — > the correct Kruzhkov's- 
type entropy solution leads to difficult open problems. We begin with: 

3.1. (2m+l)th-order NDEs for any m > 1 DO NOT generate a semigroup of 
contractions in L 1 . A first naive approach would be to try to create a standard entropy 
condition for the NDE of, say, the following form (cf. (12.41) ): 

(3.2) \u - k\ t - \ [sign(w - k)(u 2 - k 2 )] xxx < in P'(Qi) for any fcGl. 

Then Kruzhkov's-type computations with (II. ip are supposed to be performed by using his 
fundamental idea of doubling the space dimension; see a clear presentation in [23 P- 402], 
with some obvious adaptations of test functions involved. 

One should avoid doing this bearing in mind that this approach must end up with 
the contractivity property (12. 5p . which cannot be true for any PDE of order larger than 
two, since these are associated with manipulations based on the Maximum Principle for 
first-order or, at most, for second-order parabolic PDEs. This means that semigroups of 
contractions in L 1 are not available for such NDEs (11.151) with any m > 1. 

3.2. On smooth solutions and odd-order operator theory. Thus, we return to the 
Cauchy problem for the NDE (11.11) . As we have mentioned, unlike the first-order case 
(12. 6p . applying the ^approximation as in (13.11) leads to a number of principal difficult 
problems and, in the maximal generality (excluding special cases), does not give neither 
existence of a solution via the family {u e } nor uniqueness of an "e-entropy" solution, |13j . 

We will develop other concepts of solutions by different types of approximations, and 
then the concept of uniqueness will be attached to the nature of existence results. 

On local semigroup of smooth solutions. Beforehand, it is of importance that, as the simi- 
larity solutions in [HI § 3] showed, the NDE (II. ip does not admit a global in time solution 
for any bounded L 1 data. This is in striking difference with the conservation laws (12. ip . 
where such existence is guaranteed by the Maximum Principle. Therefore, we restrict our 
attention to weak solutions u(x,t) in Qi, where 

(3.3) uo(x) G C^°(M) is sufficiently small. 
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Then, as the first step of a similar construction, we have to check that for such smooth 
initial data Uq, there exists a unique local classical C x 'l solution u(x, t) of (11. ip . Recall that 
characteristic methods similar to that in (I2.8P are not available for higher-order PDEs. 
This just means that (11.11) generates a standard local semigroup in the class smooth 
functions. These results are known for non-degenerate NDEs such as (II. 6ft . and moreover 
the solutions are C°° locally in time, [21E1E2]- Actually, these smoothing results can be 
viewed in conjunction with classic methods of analytic semigroups in PDE theory; see [I] 
and references in a more recent paper [8]; see below. 

Uniqueness and continuous dependence: an illustration. Actually, in our construction, we 
will need just a local semigroup of smooth solutions that is continuous is L\ oc . The fact 
that this is generated by third-order (or other odd-order NDEs) is illustrated by the 
following easy example. Consider, for defmiteness, the NDE 

(3.4) u t = A(u) = uu xxx , u(x, 0) = u (x) G H 7 (R), 

where, without loss of generality, we take into account the principal higher-order term 
only. According to the above results, we assume that u(x,t) satisfies 

(3.5) i < u < C (C > 1) 

and is sufficiently smooth, u G L°°([0, T], H 7 {R)) and u t G L°°([0, T], H 4 (R)). See details 
on such uniqueness results in [31 § 3]. 

Thus, assuming that there exists the second smooth solution v(x,t), we subtract the 
equations and obtain for the difference w = u — v the following: 

(3.6) w t = uw xxx + v xxx w. 

We next divide (13. 6p by u > ^ > 0, multiply by w in L 2 , so, after integrating by parts, 

(o *j\ f ww t — ! _d_ f _i_ I f u t w 2 f v xxx w 2 

\°-'J J u — 2 dt J u ^~ 2 J u 2 — J u ■ 

Therefore, using the assumed regularity yields 

(3-8) U/v = /fr-|t)^ = /fr-I H r)^<C'i/^ 

where we use the fact that u xxx (-, t), v xxx (-, t) G £°°([0, T]). By Gronwall's inequality, (13.81) 
implies that w(t) = 0. As usual, this construction can be translated to the continuous 
dependence result in L 2 and hence in L\ oc . 

On degenerate NDEs. For degenerate NDEs such as (11.11) and for solutions of changing 
sign, the unique local smooth solvability is a technical result, which we do not completely 
concentrate upon, and present below some rather formal comments justifying such a local 
continuation. One of the main difficulties of this local analysis, is that (II. ip admits 
solutions with finite interfaces and free boundaries, which represent "weak shocks" with 
quite tricky (smooth enough but not C x ) behaviour. 

Thus, in addition, except the shock waves, which we are mostly interested in, the NDE 
(II. ip is degenerate at {u = 0}, so that the local existence of sufficiently smooth solution 
must include the demand of "transversality" of all the zeros (a finite number) of initial 
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data Uq(x) (or u(x,tj(5)) later on). Here the transversality of the zero at, say, x = has 
a standard meaning: 

«o(0) ^ 0. 

For instance, for key applications, we may assume that Uq(x) is anti- symmetric, so 
u(—x,t) = —u(x,t), and hence the only transversal zero is fixed at the origin x = 
only, i.e. 

(3.9) u(0,t) = 0, and u(x, t) > for x < 0. 

Then, according to regularity results for odd-order PDEs [2J El EH EH] (cf. [1 EJ [27]), the 
linearization about sufficiently smooth Uq(x) yields that the possibility of local smooth 
extension of solution is governed by the good spectral properties of the third-order linear 
operator with the principal part 

(3.10) pi = x|L f or x ~o+. 

This type of degeneracy is not sufficient to destroy good spectral properties of P3 that 
still will admit a discrete spectrum and a compact resolvent in the corresponding weighted 
space ~ Lj. for x > 0. Note that the singular point x = starts to generate a continuous 
spectrum for the operator 

(3.11) ^ = x n ^ (x>0) 

in the parameter range n > 3 only, i.e., for much stronger degeneracy than in (I3.10p . 
Indeed, then the change z = x a with a = > transforms (13. lip into the regular 
operator with the constant principal part 

(3.12) P3 = D 3 Z for z^0 + , 

for which all necessary spectral properties are obviously valid, [30]. The finite interface 
behaviour will be shown to correspond to n = 2, so it is still in the good range. Our 
conclusions here are based on the well-known fact that the linear PDE 

(3.13) u t = u xxx 

generates a smooth (analytic in a properly weighted L 2 -space) semiflow given by 

(3.14) u{x,t) = b(x — -,t) * %(•), 
where b(x, t) is the fundamental solution 

(3.15) b{x,t) = r*F(x/t%), where F = Ai(z), F" + \Fz = $ 1 fF=l. 

Thus, for the degenerate NDE ( 11.11) . the notion of "sufficiently smooth solutions" should 
also include the assumption of transversality, i.e., of local behaviour near zeros. Of course, 
this is not that essential hypothesis that has a local character, and, for instance, completely 
disappears for the related non-degenerate NDEs such as (jl.lOp . which also admits shocks 
and needs proper entropy theory (to be treated also). 
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On odd- order ordinary differential operators. In the above analysis, we need a detailed 
spectral theory of third-order (or more generally, odd-) operators such as 

(3.16) P 3 = a{z)Dl + b(z)D 2 z + c(z)D z + d{z)I, ze{-L,L) (a(z) > c> 0), 

with bounded coefficients. This theory is available in Naimark's classic book [301 Ch. 2]. 
It was shown that for regular boundary conditions (e.g., for periodic ones that are regular 
for any order and that suit us well), operators admit a discrete spectrum {A^}, where the 
eigenvalues A& are all simple for k 3> 1, and a complete in L 2 subset of eigenf unctions 
{ipk} that create a Riesz basi£|. This makes it possible to use standard eigenfunction 
expansion techniques; see necessary details and references at the end of Ch. 2 therein. 
The eigenvalues of (13.161) have the asymptotics 

(3.17) A fc ~ (±2vrA;i) 3 for all k » 1. 

In particular, this means that P3 — al for any a ^> 1 is not a sectorial operator that 
makes suspicious referring to the analogies with analytic theory [H El [27] that is natural 
for even-order parabolic flows. 

Nevertheless, recall that (I3.14p guarantees analyticity of solutions that is now associated 
with the Airy-type operator 

(3.18) B 3 = ^ 3 + l zDz + 1 1 i n L 2(R), p(z) = e a ^\ 

where a > is sufficiently small; cf. a "parabolic" version of such a spectral theory in [7J. 
It turns out that (13.181) has the real spectrum (see [121 § 9]) 

a(B) = {-|, Z = 0,1,2,...}, 

so that B — al is sectorial for a > (Ao = is simple), and this justifies the fact that 
(13.141) is an analytic flow. 

Note also that analytic smoothing effects are known for higher-order dispersive equa- 
tions with operators of principal type, [35] . This suggests to treat (13.41) by classic approach 
as in Da Prato-Grisvard [I] by linearizing about a sufficiently smooth uo = u(to), to > 0, 
by setting u(t) = u + v(t) giving the linearized equation 

(3.19) v t = A'(u )v + A(u )+g(v), t > t ; v(t ) = 0, 

where g(v) is a quadratic perturbation. Using good semigroup properties of e A/<n °)*, this 
makes it possible to study local regularity properties of the integral equation 

(3.20) v(t) = J e A '( u °)( t - s )(A(u ) + g(v(s))) ds. 

to 

It is key that the necessary smoothness of solutions demands the fast exponential decay 
of solutions v(x,t) as x — >• 00, since one needs that v(-,t) G L 2 p ; cf. [25J, where C°°- 
smoothing also needs an exponential-like decay. Equations such as (13.201) can be used to 
guarantee local existence of smooth solutions of a wide class of odd-order NDEs. 

lr This is G.M. Kessel'man's (1964) and V.P. Mikhailov's (1962) result. 
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Thus, we state the following conclusion to be used later on: 

any sufficiently smooth solution u(x, t) of (13.4)) . (13. 5p at t = t , 
can be uniquely extended to some interval t E (t ,t + 5), 5 > 0. 

3.3. Global solutions by Galerkin method. Here we demonstrate the application 
of another classic approach to nonlinear problems that, suddenly, in the present case of 
unclear entropy nature of solutions of NDEs and the open uniqueness problem, gives a 
partial answer to both. We mean the Galerkin method that was the most widely used 
approach for constructing weak solutions via finite-dimensional approximations; see Lions 
[26J with many applications therein. 

Thus, by this classic theory of nonlinear problems, under the assumption (13. 3p and 
others, if necessary, let us perform a standard construction of a compactly supported (for 
simplicity) solution by Galerkin method using the basis {ipk} of eigenfunctions of the 
regular linear operator P 2 = D 2 < with the Dirichlet boundary conditions, 

(3.22) ip" = Xkip, ip — at x = ±L ==>- ~ —k 2 . and u = at x = L. 

As an alternative, it is curious that, for our purposes, possible (and more convenient for 
some reasons) to use the eigenfunction set of the operator P 4 = — < again with the 
Dirichlet conditions 

ip — if>x = at x = ±L. 
Special Galerkin bases associated with higher-order operators Pq = D® < are also may 
be convenient; see applications to third-order linear dispersion equations in [21]. 

In all these self-adjoint cases, the eigenfunctions form a complete and closed set in L 2 ; 
see classic theory of ordinary differential operators in Naimark [3U1 P- 89]. 

On the other hand, looking more natural choice of the third-order operator P3 = D\ 
for Galerkin approximation of (II -ip will cause a difficult problem, since for the third-order 
PDE with the principal operator as in (jl.ip . 

(3.23) u t = a(x, t)u xxx + ... (ohm) 

with a > 0, proper setting for the IBV problem includes the Dirichlet conditions (see 
Faminskii [TT] for details and a survey) 

(3.24) u = u x = at x — —L and u = at x = L. 

For a < 0, the boundary conditions must be swapped, so that the proper setting of 
the problem depends on the unknown sign of solutions. Here, the fact that P3 = D% 
is not self-adjoint is not essential since, relative to adjoint basis {if>l}, the closure and 
completeness of the bi-orthonormal generalized eigenfunction sets remain valid. 

Actually, the choice of linear operators P 2 = D%, P4 = —D*, or others, is not of 
principal importance if we are looking for compactly supported solutions 

(3.25) u E C^°((—L, L) x [0, 1]). 

It should be noted that the control of finite propagation property in (11.11) is difficult and 
is an essential part of our further analysis. For instance, we also can fix periodic boundary 
conditions that are always regular, [301 Ch. 2] (it is curious that (I3.24p are not). 
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Thus, we construct a sequence {u m } of approximating Galerkin solutions of (jl.ip . (11.21) 
in the form of finite sums 

m 

(3.26) u m (x,t) = C k {t)4j k (x), 

k=l 

where {Cj} solve the quadratic dynamical systems 

(3.27) Cj = Yj CkQJkij, where J kij = {if>ki>'i, V>" ) = ^ji^k^i, 

(W) 

For the conservation law ( 12. lft . the DS takes the same form as in ( 13. 2 7ft . with the only 
difference that 

(3.28) Jkij = -(^'i,^). 
The identity (11.91) for u m takes the form 

(3-29) Ep^C'IW = Co m = E P ^C' fc 2 (0), t>0. 

(fc) 1 1 (fe) 1 1 

This guarantees global existence of the solutions u m (x,t) showing that 

(3.30) Ck(t) do not blow-up and exist for all t > 0. 

Since ipk are given by sin(Afcx) or cos(Afcx), a lot of coefficients J k ij vanish. For instance, 
if uq(x) is odd, we take all the sin-functions, 

^(x) = ^sin(^f), with A fc = -££, k — 1,2, ... . 

The non-zero coefficients J k ij occur iff k = j, I = 2j, where (I3.27P becomes simpler, 

(3.31) C^^-CjCy, j>l. 

It is curious that (13.311) yields the following feature of a "maximum principle" : 

(3.32) signC^t) = sign (7,(0), j > 1. 

Other a priori estimates are obtained by multiplying (II. ip in L 2 by u and u xx yielding 
the identities 

(3-33) \iJv? = -\S{u x )\ \f t !{u x f = -\ju x {u xx f. 

Then some interpolations of various terms in the identities ( 13.331) are necessary. 

Thus, the sequence of "regularized" solutions (Galerkin approximations) {u m (x,t)} is 
globally defined, and 

(3.34) {u m } is uniformly bounded in L°°([0, 1]; H' 1 ). 

Therefore, along a subsequence, {u m } converges to u weakly-* in L°°([0, 1]; H~ l ), and, in 
addition, strongly in if _1 ([0, 1]; H~ 2 ), in view of compact embedding. This gives a weak 
solution. As usual, the better regularity comes from the special choice of Galerkin's basis 
employed. We do not stress attention to this (bearing in mind local C^-smoothing for 
non-degenerate NDEs). See [23] for rather exotic Galerkin bases applied to KdV type 
equations. Recall that, globally, smoothing is not available, since this construction is 
specially oriented to include shocks of SL-type. 
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Remark 1. Obviously, the estimate (I3.29P does not and cannot prevent gradient catas- 
trophe, which means that 

(3.35) \\ux(t)\\ 2 2 = E \ X k\Cl(t) -> +00 as t -> T~ < 1. 

Notice that for (11.11) there is an opportunity to create blow-up of the solutions u(-,t) itself 
(possibly together with (13.351) ). where 

(3.36) \\u(t)\\l = J2 C k(t) -»• +°° as t^T"<l. 

This does not happen if a finite shock appears via the self-similar patterns such as [18] 

(3.37) u-(x,t)=g(z), z = x/(-t)*, {gg')" = \g'z, /(too) = ±1. 
Indeed, by the first identity in ( 13.331) . there appears an integrable singularity, 

(3.38) \± /(«.)» ( / 7 (ffO»<keLi((-l,0)), 

-1 -1 -(-t)- 1 / 3 

so that ||w(0 - )||| remains finite. Here in (13.381) one needs to use the asymptotics of the 
Airy function [T8"j § 3], so that the integral therein diverges but its rate, 

(-t)-V3 

/ (g'fdz ~0((-f)-*), 

-(-t)-V8 

is sufficient for the integrability. 

Remark 2. Using the dynamical system (13.271) instead of the NDE (II. ip suggests to 
develop a formal calculus of the corresponding sequences, where, on identification, 

(3.39) u = ECW>* u = {C k } 
belongs to the little Hilbert space hp 1 with the metric 

(3-40) INI^E^ 2 . 

Then (I3.29P guarantees that 

(3.41) u{t) E hp 1 for all t > 0, 

meaning global solvability. Moreover, the embedding hp 2 C hp 1 is compact since t^j ~ p 
[28] (for Zip 2 , the metric contains in (13.401) ). so that we can use the same Galerkin ap- 
proximation method to construct suitable solutions. In this space, the blow-up formation 
of shocks means (I3.35p . 

Remark 3. Writing the iV-dimensional NDE (11.161) for compactly supported Uq as 
(3-42) (-A)-V = -^ 2 

with the standard definition of the linear operator (— A) -1 in L 2 (Q), Q is sufficiently large, 
and multiplying (I3.42p by u yields the same conservation identity (II. 9p . Some concepts 
developed above can be also adapted to the equations in Mr, though shock wave formation 
phenomena become more involved and are in general unknown. 
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3.4. ^-entropy solutions (a test) for the NDE. Assuming that the local smooth 
solvability problem above is well-posed, we now present the corresponding definition that 
will be applied to particular weak solutions. Recall that the topology of convergence, L\ oc 
at present, for (11.11) was justified by a similarity analysis presented in [TBI Prop. 3.2]. For 
other NDEs, the topology may be different that can be a difficult problem. 

Definition [3jl. A weak solution u(x,t) of the Cauchy problem (11.11) . fl X . 2 1) is called 
5-entropy if there exists a sequence of its smooth 5- deformations {u Sk , k = 1, 2, ...}, where 
5k — > 0, which converges in L\ oc to u as k — ► oo. 

Note that this is slightly weaker (but equivalent) to the condition of 5-stability. The 
construction of global 5-deformation of u is performed along the lines of (i)-(k) in Section 
12.21 The only difference is that local <5-deformations can lead to complete blow-up for 
the NDE ( 11. ip . as explained in [181 § 4.2]. To avoid this, one needs either to impose 
the condition (13.31) or specially assume that complete blow-up cannot occur under slight 
deformation of the data, or while performing its 5-deformation with any sufficiently small 
5 > 0. We call such solutions 5 -extensible (the definition assumes that u is (^-extensible). 

On 5-entropy test and uniqueness. First of all, we again note that any uniqueness (and 
entropy) results for such NDEs are not acievable in principle, [15]. Therefore, we use the 
above results as a basis of the so-called " 5- Entropy Test" for testing shock an rarefaction 
waves; see first applications below. 

S-entropy solutions: motivation of the term. Let us explain why solutions are called 5- 
entropy, while we do not use any evolution integro-differential inequality such as (12.41) . It 
turns out that the NDE (II. ip itself contains the right evolution choice of the admitted 
type shocks in the class of smooth solutions (precisely this makes sense of Definitional). 

For instance, as a rough explanation, assume that at x = the shock S+ is going 
to appear at t = l - from a smooth solution u(x,t) such that u(x,l~) remains smooth 
everywhere except x = 0; e.g., for simplicity, we assume that 

(3.43) u(x, 1) ~ S+(x) in a neighbourhood x G (—5,6), 

together with necessary derivatives u x and u xx that are assumed to be small at x — ±5. 
Here 5 > is also a small constant, so our illustration is of local nature. Multiplying (11.11) 
by u and integrating over (—6, 5) for t ~ 1~ yields the following main terms: 

(3.44) l± j{u 2 - 1) = \j(u x f + ... > (or < for at t = l), 

since u x must be essentially positive on profiles u(x,t) that smoothly approximate S + (x). 
One can see that ( 13. 44ft evolutionary prohibits stabilization to S + (x) as t — > 1~ , when 
v? — ► 1 in L\ oc . More rigorously [T31 § 7.2], the same negative result is established using 
the weaker topology of H~ l , where multiplication applies to the non-local equation (jl.8p . 

Similarly, we arrive at no contradiction while using (I3.44p to describe stabilization to 
S-{x), since then u x is essentially negative. In fact, (13.441) reflects a finite-time formation 
of the singular shock S- (the gradient catastrophe) for the NDE (11.11) that was described 
in [TU § 3] in greater detail. 
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Thus, using smooth deformations guarantees (via smoothness, that is important) the 
preservation of the natural local entropies such as inequalities like fl3.44j) and the opposite 
one for S-, so we call the constructed solutions 5-entropy. 
First easy application of 5 -entropy test. As a first application, we have: 

Proposition 3.1. Shocks S_(x) and H(—x) are 5 -entropy. 

The result follows from the properties of similarity solutions (13. 3 7ft . which, by shifting 
the blow-up time T f— > T + 5, can be used as their local smooth 5-deformations at any 
point t E [0, 1). For H(—x), we will need an extra approximation of similarity profile g(z) 
with finite interface at some z = z , at which it is not C 3 , by sufficiently smooth profiles. 

Let us use the negation in the following form: 

Definition [3J.2. A weak solution u(x,t) of the Cauchy problem (II. ip . (11.21) is not 5- 
entropy if it is not 5-stable. 

Proposition 3.2. Shocks S + (x) and H(x) are not 5-entropy. 

Indeed, taking initial data S+(x) and constructing its smooth 5-deformation via the 
self-similar solution [18, § 3.4] 

(3.45) u+{x,t) = g[z), z = x/t*, {gg')" = -\g' 'z, /(Too) = =f1. 

Performing time-shifting t i— ► t+5, we obtain the global 5-deformation {u 5 = u + (x,t + 5)} 
which goes away from S + . 

Thus, we have shown that, at least, the idea of 5-deformations allows us to distinguish 
basic 5-entropy and non-entropy shocks without any use of mathematical manipulations 
associated with standard entropy inequalities, which are illusive for higher-order NDEs 
(and nonexistent in principle [T5]). 



4. COMPACTONS ARE 5-ENTROPY SOLUTIONS 

Without loss of generality, we treat this question for a particular NDE. Namely, consider 
the following K(2,2) equation: 

(4.1) Ut = (uu x ) xx + 4uu x in Kx M. + . 
Its compacton solution has the explicit form [33J, 

(4.2) «.(,,*)=/„(, + *), where f M = { 2c f «> 

This is an example of a compactly supported weak solution of equation 04. ip . One can see 
that at the interface points y = ±tc, the profile f c (y) is just C*' 1 , i.e., the first derivative 
f' c (y) is Lipschitz. Therefore, it is not a classical C x 'l solution of the PDE and has weak 
singularities at y — ±7r, so one needs to check whether it is an entropy solution. In 
addition, the "flux" (//')' is continuous at those points, though this does prove nothing. 
We now use the concept of 5-entropy solutions from Sections [2] and [3j 



Proposition 4.1. The compacton {4-2) is a 5-entropy solution of the NDE (4-1)- 
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Proof. We are going to show that there exists smooth 5-deformations of u c for arbitrarily 
small 5 > 0. The general TW solutions as in (14. 2p with A = —3 yields the ODE 

(4.3) 3 f = (//')" + 2(/ 2 )' => 3/ = (//')' + 2f 2 + C 5 , 
where we chose the constant of integration to be 

(4.4) C s = 35 - 25 2 > 0. 

One can see on the phase plane in the variables {f 2 , (f 2 )'} that the ODE (14.31) . (I4.4p has 
a strictly positive and hence analytic solution fg satisfying 

(4.5) fs(y) - ^ + as V —* =t°° exponentially fast, and/5 — > f c as 5^0 + 

uniformly in R. According to Definition [3)1, (14.51) implies that u c is an entropy solution 
of (14.11) . as well as f c is G-admissible for the third-order ODE in (14.31) . □ 

It is crucial that Proposition 14.11 justifies that the K(2,2) equation (14. ip . the NDE 
(II. ip . and many others with similar principle degenerate third-order operators possess 
finite propagation of interfaces for entropy solutions. 

It is worth recalling again that, regardless the existence of such nice smooth compactons 
(I4.2p . the generic behaviour for the RH equation (14.11) . for other data, includes formation 
of shocks in finite time, with the local similarity mechanism as in [TBI § 3.1]. 



5. On extensions to other related NDEs 

5.1. Shocks for the non-degenerate NDE. We begin with the simpler model (11.101) 
that appeared in Section [3T21 while we discussed the possibility of extensions of sufficiently 
smooth solutions for defining (^-deformations. Indeed, for (11.101) . this is much easier. On 
the other hand, obviously, as an NDE, this admits shocks via standard similarity solutions 

(5.1) u.(x,t)=g(z), z = x/(-t)\ => {{l + g 2 )gi)"=\g'z. 

This ODE is studied as usual. Figure [T](a) shows a few similarity profiles satisfying 

(5.2) g(z) ~ z~4e~ a ° z3/2 — > as z — > +00, where a = ^4g, 

that create as t — > 0~ the shocks ~ H(—x). By dotted lines, we indicate there other 
profiles g{z), for which g(+oo) 7^ 0. For the sake of comparison with compactons, in 
Figure \Bjo), we present the soliton of the related NDE 

(5.3) u t = ((1 + u 2 )u x ) xx + (1 + u 2 )u x , where 

(5.4) u s (x,t) = Uy), y = x-Xt => -\f ={(l + f 2 )f)" + (l + f 2 )f. 
The soliton profiles have now exponential decay for A < — 1, 

fM ~ z- aoU -> as |yH+oo, a = v / \l + X\. 
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Similarity shock profiles of type H(-x) for u t =({1+u 2 )u x ) xx TWsoliton, u=((1+u 2 )u x ) xx +(1+u 2 )u x ; X=-2 




-100 -80 -60 -40 -20 10 15 20 25 30 35 40 



z y 

(a) shock profiles (b) soliton 

FIGURE 1. Similarity profiles for the NDE (jl.lOp : shock profiles satisfying the 
ODE in (pUj) (a), and TW soliton satisfying ([53]) with A = -1. 

5.2. (5-entropy approach to the NDE— (2,1). For the non-fully divergent PDE (11.121) 

we also apply the (5-entropy to prove existence and uniqueness via suitable approximations. 

On Galerkin method . Constructing Galerkin approximations, we face a new technical 
difficulty in passing to the limit since a uniform estimate such as (13.291) is not available 
for solutions [18, § 4] 

(5.5) u a {x,t) = (-t) a g{z), z = x/(-t) f) , p = ijp (a G R). 

Nevertheless, we can establish some extra estimates by using the corresponding DS (I3.27p . 
where Jkij = —Xi{ipk4'i,'4 , 'j)- E.g., for odd data, the simpler system similar to (13. 3 1 j) . 

Cj = jofCjCv, 1 < j < m, 

implies that, for m even, 

(5.6) CL = 7o (f ) 3 C f C m and Cj = 0, j > f C f (i) = C f (0)e^(f ) 3c -(°)*. 
Therefore, assuming that 

(5.7) M0 GC 3 |C m (0)|<^r, m>l (c>0), 

we obtain from (15.61) a uniform bound on the Galerkin coefficients {Cj}, and hence a local 
weak solution. 

Shocks and compactons exist. On the other hand, regardless its non-full divergence and 
nonexistence of any obvious conservation laws, the NDE ( 11.121) allows a similar treatment 
of shocks and rarefaction wave as for ( 11.11) . For instance, formation of finite shocks for 
( 11.121) is described by the same self-similarity as ( 13.371) . with the ODE, 

(5.8) u-(x,t) = g(z), z = x/{-tf* => {gg")' = \g'z, g( T oo) = ±1. 

Existence and uniqueness for (15. 8p is proved similar to [TH Prop. 3.1]. In Figure [2(a), we 
show a few similarity profiles that create as t — > CT the shocks. The profile for S-(x) has 
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m=1 , NDE-(2,1 ): basic shock similarity profile for data of S_-type 




m=1, NDE-(2,1): the shock similarity profile for H(-x), z n =1.297... 

g'T=g'2/3 I 



10 20 30 40 




(a) for S- (x) (b) for H (-x) 

FIGURE 2. The ODE (|5.8p : the shock similarity profiles including the unique 
solution (boldface line) for data S-(x) (a), and for H(—x) with finite right-hand 
interface (b). 

m=1 , NDE-(2,1): saw-type shock similarity profile for a =-0.2384... 




-40 -30 -20 -10 



10 20 30 40 



FIGURE 3. The saw-type similarity solution of the ODE in ([579]) for a c = -0.2384... . 

the derivative at the origin 

g '(0) = -0.702... (instead of g'(0) = -0.51... for the NDE-3 (O]) .) 

In (b) explaining formation of H(—x), the right-hand interface is situated at zq = 1.297. 
As another known key feature, Figure [3] shows the saw-type profile for the ODE 



l + a 



(5.9) 



u-(x,t) = (-t) a g(z), z = x/(-t)~ 



(99")' = 1 -^9'z-ag 



where a. 



-0.2384. 
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Similarity shock profile dissappears as n tends to 2 ; g'(0)=-0.51 




-40 -35 -30 -25 -20 -15 



FIGURE 4. Shock similarity profiles of the ODE f|5.13|> for various n £ (0,2); 
g(0) = g»(0) = 0. 



The compacton equation associated with (II ,12ft takes the form 

Ut = (uu xx ) x + 2uu x 
and admits the TW solution with the same f c as in (14,21) . but now for A = —1, 



(5.10) 



u c (x,t) = f c (x + t), f c (y) 



2 cos 2 (|) for \y\ < ir, 



for \y\ > 7T. 

As for (11,11) . it is 5-entropy and G-admissible; Proposition 14.11 is proved similarly. 

5.3. Shock similarity profiles for Harry Dym-type equations. Consider the NDE 
(5.11) u t = \u\ n ~ 1 uu xxx (n>0), 

which for n = 3 becomes the quasilinear Harry Dym equation 



(5.12) 



Ut — u u a 



which also belongs to the NDE family and is an exotic integrable soliton equation; see 
§ 4.7] for survey and references therein. It admits the same formation of shocks S_(x) by 
the similarity solutions given in (I5.8P with the ODE 



(5.13) 



\9\ n - 1 99'" = i9 f z. 



Figure H] shows that such similarity profiles exist for n G (0, 2) and vanish as n — ■> 2 
(proof is easy), so that for n = 3 (the Harry Dym case) such shocks are not available. 
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Shock similarity profiles for fully nonlinear NDE |u t |' u t =(uu x ) xx 




FIGURE 5. Shock similarity profiles of the ODE in (j5.14jl for various 7 G (—1, 4]; 
g(0) = g"(0) = 0, g'(0) = C = -0.51. 



5.4. Shocks for fully nonlinear NDE. For the NDE (I1.13p . the basic blow-up similar- 
ity solutions are slightly different, 



(5.14) u_(x,t) = g(z), z 



(99')" = P 1+ ~ f \g'z\^g'z. 



Mathematics of such ODEs is not much different than that for (I3.37[) . In Figure [5], we 
show how the shock similarity profiles g(z) depend on 7 > —1. All these profiles satisfy 
the anti-symmetry conditions at the origin, 

(5.15) g(0) = g"(0) = 0, 
and the following expansion holds: 

(5.16) g{z) = Cz + f3 1+ ~> C1 



zpz 3 + ... . 



(2+7) (3+7) (4+7) 

Note that the linearization about the constant equilibrium C. 
yields a nonlinear ODE, 



1 as z 



-00, again 



(5.17) 
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f3 1+1 \g'\ 1 g'\z\~<z + ... 



which deserves further study. Figure [5] shows that the solutions remain equally oscillatory 
for all 7 > —1, i.e., this is not a manifestation of the oscillatory character of the linear Airy 
function that occurs at a single simplest value 7 = only. Thus, all ODEs ( 15. 14ft with 
7 > — 1 contain a strong nonlinear mechanism of oscillations about constant equilibria. 

5.5. Shock similarity profiles for cubic NDEs. Analogously, in a similarity fashion, 
the shock formation is studies for the cubic fully divergent NDE (II. lip . The formation of 
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shocks H(—x) is described by the similarity solutions (15. 8p . where 
(5.18) (9 2 97 = l9'z, 

which admits a similar rigorous study. Figure [6] shows similarity profiles with the finite 
interface at z = Zq > with the expansion as z — > zq 

(5-19) g(z) = Jf(z - z) + - i^L. (z - z)\ + ... , 



10 

for which the flux (g 2 g')' = | {g 3 )" is continuous at z — z Q , so these are weak solutions. 
The flux is not zero for a more singular expansion such as 

g{z) =C(z -z)l + ^{z Q -z)l + ... {z<z , C>0). 

Similar to [HI § 3.1], such blow-up similarity solutions describe the generic formation of 
shock waves of the type ~ H(—x) for (11.111) . These solutions are entropy, which is proved 
by regular analytic approximations of the ODE as in Section HI 

By dashed lines in Figure El we denote other profiles, for which C + = g(+oo) > 0, 
so that the corresponding blow-up similarity solutions (I3.37P lead to more general shocks 
with different values C± as z — > ±oo (with C + > 0). Then, as z — > +oo, g(z) approaches 
C + exponentially fast, 

g(z) = C + + 0(e~ a ° z3/2 ), where a = 

Thus, the above solutions with the behaviour (I5.19P close to interfaces show finite 
propagation for the NDE (11.111) . There are also TWs with finite interfaces given by 

u (x,t) = f(x+t) =>■ f = (f 2 fy 

that are entropy and are approximated by the analytic family {/,5 > 5 > 0} satisfying 

f-5 = (f 2 f')' (5>0). 

For instance, the following TW with the interface at y — is 5-entropy: 

^ {-y) for y < 0, 
for y > 0. 



f(y) 



Other discontinuous TWs may not admit smooth approximations via similar TWs. 

The boldface line in Figure [6] indicates the profile that leads to H(—x) as t —>■ 0~. Here 
the shock H(—x) is not a weak solution of the NDE (15.181) . Recall that it is a 5-entropy 
solution, i.e., there exists a converging sequence of its smooth (^-deformations. 

More advanced shock patterns are created by similarity solutions (15.51) . with 

(5.20) /3 = ±f* and {g 2 g')" = ±±p g'z - ag. 
The interface expansion (I5.19P changes into 

(5.21) 9W = ^^ (zo _ z)+ _i 7 _i= 
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z - z) 2 , + 



Shock formation similarity profiles for u t =(u 2 u x ). 




-60 -50 -40 -30 -20 -10 10 

Z 

FIGURE 6. Shock similarity profiles satisfying (j5.18j) . (15.19P for various Zq > 0; 
the boldface profile leading to H(—x) has zq = 1.20...; dotted lines denote shock 
profiles with g(+oo) > 0. 




FIGURE 7. Shock similarity profiles of the ODE in (|5.20p satisfying expansion 
(|5.2ip for zq = 1.2 for various positive and negative a. 

Figure [7] shows typical solutions of the ODE in (15.201) for a > and a < 0. The most 
interesting "saw-type" profiles occurs at 



a c w -0.0715. 
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Shock formation similarity profiles for u ( =(u u xx ) x 




FIGURE 8. Shock similarity profiles satisfying the ODE flgjUP - 
On non- divergent cubic equation. Consider briefly the cubic NDE-(2,1), 



(5.22) 



u t 



{u 2 u. 



which is similar, though it does not admit finite propagation at the degeneracy level 
{u = 0}. This is seen by using TWs 



(5.23) 



so, on integration twice, 

Ci 



u(x,t) = f(x+t) => r = (/ 2 /t, 

f" = i + ^ =► Hf') 2 = ^\f\-f + c 2 . 

Setting C\ = by assuming continuity of flux: f 2 f" = at / = 0, yields the ODE 

i(/') 2 = ln|/|+C 2 

that does not allow any connection with the singular level {/ = 0}. 

The shock similarity profiles for ( 15.221) exhibit the same form (15.81) and the ODE is 



(5.24) 



(9 2 9 



2„ll\l _ 1 _/_ 

- 3 y z. 



Typical strictly positive profiles with g(— oo) = C_ > C + = g(+oo) > are shown in 
Figure [U so these describe blow-up formation of more general entropy shocks. 
More general blow-up similarity patterns (15.51) for (15.221) yields the ODE 



(5.25) 



(g 2 g"Y = l -^g'z- ag, 



which exhibits properties that are similar to (I5.20p . In Figure [9^a), we show typical 
solutions of (I5.25P for a = — These profiles are strictly positive with 



(5.26) 



g(z) ~ \z\ !+ 2q as z — >• ±oo (a G (— |, 0)). 
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Shock similarity profiles for u =(u u ) ;a=-0.1 



"Saw-type" shock similarity profile; a =-0.1 2559.. 
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(b) "Saw" for a c = -0.12559 



FIGURE 9. The ODE (f5T25"j) : the shock similarity profiles for a = (a), and 



a "saw" profile for a = a c = —0.12559... (b). 

In the bottom right-hand corner of (a), we present a number of "steep" solutions that 
quickly vanish (according to ( 15.231) with C\ < 0). These show that the asymptotics (15.261) 
is unstable in the direction of shooting from z = +oo. 

In (b), we present a special profile that plays a role of the "saw- type" solution for 



a 



a,. 



-0.12559.... 



u c (x,t) = f c (x + 2t), f c (y) 



This is the best "saw" we can get numerically, though it is seen that there exists the first 
vanishing point while other "teeth" still stay away from zero. Anyway, we have checked 
that positive shock profiles cannot be extended to a < a c , so this is definitely a critical 
value of parameter. 

Related compactons. Consider the following compacton equation (q.v. (11.111) ): 

u t = {u 2 u x ) xx + 9u 2 u x . 
The explicit compacton solution is now easier, 

cosy for \y\ < ~, 
for \y\ > §. 

Regardless the fact that it is not C 1 at the interface, this solution is ^-entropy (note that 
(15.191) exhibits the same regularity). The proof uses regular approximations as in (14. 3p . 

5.6. An analogy with parabolic problems. In a natural sense, an analogy of the 
difference between the NDE-3 (NDE-(0,3)) (II. ip and (11.121) can be observed in nonlinear 
parabolic theory. Namely, the fully divergent fourth-order diffusion equation (the DE-4, 
or DE-(0,4)), 

(5.27) u t = -{\u\u) xxxx in Q x = (-L, L) x (0, 1) 

(recall that the nonlinearity \u\u keeps the parabolicity on solutions of changing sign), by 
classic parabolic theory [261 Ch. 2], admits a unique weak solution of the Cauchy-Dirichlet 
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problem with data u such that (u ) 2 G H 2 . Multiplying (I5.27P by in L 2 {Qi) and 

integrating by parts yields the following a priori estimates of such weak solutions: 



\u\u G L°°(0, T; H ) and {y\u\ u) t G L (Qi). 

Uniqueness follows from the monotonicity of the operator — (\u\u) xxxx in H~ 2 : for two 
weak solutions u and v , 



= — f (\u\u — \v\v)(u — v) < 0, 

so that (I5.28P guarantees continuous dependence of solutions on initial data. 
On the other hand, the fourth-order thin film equation (TFE-4) 

(5.29) u t = -(\u\u xxx ) x , 

which has the distribution of the derivatives (3,1), does not admit such a simple treatment 
of continuous dependence and uniqueness as via (I5.28p . The Cauchy problem for the non- 
fully divergent TFE-4 (15.291) needs special approximation approaches, [9]. 

For non-fully divergent operators such as in (I1.12p or fifth-order ones of the types (2,3), 
(3,2), (4,1), in the NDEs (see [T4]) 

Ut {nu xx } xxx , Ut {nu X xx)xxi ^t {uu xxxx } x ^ 

we face a difficulty that is similar to that for the TFE (15.291) . In both cases, the 5- 
approximation concepts will play a role, quite similarly to the higher-order parabolic 
TFEs-6 such as (see [10] and references therein) 

Ut ( \u \u XXXX x) xi Ut (|w|ti^^ :r:r ) a ; a ;, Ut (|w|'U :2 ;a;a;)a;a;a:) etc. 

6. On related higher-order in time NDEs 

It is principal for PDE theory to justify that the ideas of similarity shock wave formation 
remain valid for other NDEs that are higher-order in time. We claim that the concept of 
smooth 5- deformations can be developed for such quasilinear degenerate PDEs. Let us 
present a few comments in these directions. 

6.1. Second-order in time NDE. As in [TBI § 1-2], we begin with the simple observa- 
tion: S±(x) are stationary weak solutions of the second-order in time NDE 

(6.1) u tt = (uu x ) xx . 

To distinguish the entropy one, as usual, we introduce the similarity solutions 

2 

(6.2) m_(x, t) = g(z), z = x/(—t)z, where 

(6.3) {ggr = f g > z + l g » z 2 in Rj /( Too ) = ±i. 

The study of this ODE is similar to that in [TH § 3] , so we present the existence result for 
the shock S-(x) in Figure [TOT The dotted lines show nonexistence of similarity profiles for 
S+(x) (cf. a proof below). The boldface profile is unique and satisfies the anti-symmetry 

25 



Shock formation similarity profiles for u tt =(uu x ) xx , g'(0)=-0.57. 
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FIGURE 10. The shock similarity profile as the unique solution of the problem (|6.3 



conditions at the origin (I5.15p . We see that profiles g(z) are now non-oscillatory about 
±1 and the convergence to these constant equilibria is exponentially fast, 

g(z) = ±1 + 0(e~27 N 3 ) as z _> zpoo. 

This reflects the fact that the fundamental solutions of the corresponding linear PDE 

(6.4) u tt = u xxx 

is not oscillatory as x — > oo. Obviously, the blow-up similarity solution (16.21) generates in 
the limit t — > 0~ the shock S-(x), i.e., 

(6.5) u-(x, t) — >• S-(x) as t — > CT. 

In Figure [111 we show various non-symmetric shock similarity profiles with different limits 
as z — >• ±oo. 

Incidentally, S+(x) cannot be obtained in such a limit, since the ODE (16.31) does not 
admit suitable similarity profiles g. This can be seen from the identity obtained by 
multiplying the ODE (16. 3j) by g' and integrating over (0, oo) with conditions (15.151) . 



-U9W 3 -Jg(9") 2 = U(9Tz>o. 



3 





-oo, since the 



Therefore, for g'(0) > 0, there is no positive solution g(z) — > +1 as z - 
left-hand side is then strictly negative. 

As an important difference with the previously studied PDEs, the NDE (16.11) is sym- 
metric under the time-reflection, 



(6.6) 



t ^ -t 
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Shock formation similarity profiles for u ={uu J 



Non-symmetric shock similarity profiles for u =(uu ) 
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FIGURE 11. Non-symmetric shock profiles satisfying the ODE in (|6.3 



so that the stationary shock S-(x) that appears as t — ► 0~ according to the similarity law 
(I6.2p (according to "centre/stable manifold" behaviour in [T8t § 6]) will next disappear in 
the same smooth similarity manner (16. 2p . where (— t) is replaced by t. As a next step, the 
concept of smooth (^-deformations should be applied to (16. ip to produce a unique solution 
of the Cauchy problem, but this demands extra more technical study. 

6.2. Third-order in time NDE. Consider the third-order in time NDE 
(6.7) uut = (uu x ) xx , or 



(6.* 



being a first-order system with the characteristic equation A 3 = u, with one real and two 
complex eigenvalues for u 7^ 0, so it not hyperbolic. 

Quite analogously, S±(x) are stationary weak solutions of (16. 7p for which the basic 
(with a = 0) similarity solutions are 




(6.9) 



U-(x, t) — g(z), z — x/(—t), where 



(6.10) {gg')" = (z 3 g')" = 6g'z + 6g"z 2 + g"' 'z 3 in E, /(=Foo) = ±1. 



Integrating (16.101) twice yields 

gg' = z 3 g' + Az + B, with constants A, B el, 
so that the necessary similarity profile g(z) solves the first-order ODE 
(6.11) S = ^ where ^4 = G/(0)) 2 > 0. 



By the phase-plane analysis of (16. lip , we easily get the following: 
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Proposition 6.1. The problem (I 6. 1 (ft) admits a unique solution g(z) satisfying the anti- 
symmetry conditions (15. that is positive for z < 0, monotone decreasing, and is real 
analytic. 

Such basic anti-symmetric similarity profiles are shown in Figure [121 These satisfy the 
expansion near the origin, as z — > 0, 

(6.12) g(z) = £ c k z 2k+1 = Cz + \z 3 + l 2 ±z 5 + ... (C = g'(0) < 0). 

(fe>0) 

Substituting the expansion in (16.121) into (16. lip yields 

g\g-z 3 )=Az =s> £ (2k + l)c k ( Cj - 6 n )z 2 ^ +1 = Az. 

(k,j>0) 

The corresponding algebraic system for the expansion coefficients {c^} is uniquely solved 
giving the unique analytic solution. The boldface profile g(z) in Figure [12] (by (16.91) it 
gives S-(x) as t — > 0~) is non-oscillatory about ±1 with the algebraic convergence 

g(z) = ±1 + j + ... as z — > =foo. 

Again, the fundamental solutions of the corresponding linear PDE 

(6.13) u m = u xxx 

is not oscillatory as x — > ±oo. The linear PDE (16.131) exhibits some finite propagation 
features with the corresponding test consisting of checking the TWs, 

u(x, t) = f{x - Xt) -A 3 /'" = /"', i.e., A = -1, 

where the profile f(y) disappears from. This is similar to a few other well-known canonical 
equations of mathematical physics such as 

Ut = u x (dispersion, A = —1) and u u = u xx (wave equation, A = ±1). 

Any finite propagation is not true for (16. 4p . The blow-up solution (16.91) gives in the limit 
t — > _ the shock S-(x), and (I6.5P holds. In Figure [T2l we also show the results of shooting 
with g'(0) > giving unbounded profiles g(z) ~ z 3 as z — >• ±oo. As usual, this means 
nonexistence of similarity blow-up profiles corresponding to 5 + -type shocks. 

A key difference with the previous problems is that the original ODE (16.101) written as 

(6.14) (g - z 3 )g"' = 6g'z + 6g"z 2 - 3g'g" 

has, instead of {g = 0}, another singular line (a kind of nonlinear u light cone") 

(6.15) L : g(z) = z 3 . 

Then, formally, the existence of global solutions of (I6.14p depends on the possibility 
of a continuous transition through it. The simpler integrated form (16.111) shows that 
typical solutions do not cross L (except at the analytic point z = 0), so that "weak 
discontinuities" do no occur. 

Since (16.71) has the same symmetry 



(6.16) 



u i — > —u, 
t ^ -t, 
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FIGURE 12. The shock similarity profile as the unique solution of the problem (|6.10p . 

as (11.11) . similarity solutions (16. 9p with — t \— ► t and g(z) i— > g(— z) also give the rarefaction 
waves for S + (x), as well as other types of collapse of initial non-entropy discontinuities. 

Using the known asymptotic properties of blow-up similarity solutions (16.91) and those 
global with —t i— > t, for convenience, we formulate the following 

Proposition 6.2. The Cauchy problem for the equation l\6.7\) admits: 



(i) an analytic solution m_(x, t) in Qj- = I x (0, T) that converges as t — > T~ to the shock 
S-(x) in L\ oc and a.e., and 

(ii) for non-analytic singular initial data as t — > + given by 

(6.17) u(x,t) ^ S + (x), Ut (x,t)^± ^ Utt (x,t)^ 3 -^ 

with uniform convergence as t —>■ on any compact subset from R \ {0} (and in Ll oc for 
u(x, 0)), there exists an analytic solution in M x R + . 

Analytic 5 -deformations by Cauchy-Kovalevskaya theorem. Eventually, we start to deal 
with the third-order in time NDE (16. 7p that turns out to be in the normal form, so it 
obeys the Cauchy-Kovalevskaya (C-K) theorem p. 387]. Hence, for any analytic initial 
data vJx,0), u t (x,0), and u tt (x,0), there exists a unique local in time analytic solution 
u(x,tp. Thus, (16.71) generates a local semigroup of analytic solutions, and this makes it 
easier to deal with smooth 5- deformations that always can be chosen to be analytic. On 
the other hand, such nonlinear PDEs can admit other (say, weak) solutions that are not 
analytic. Actually, Proposition 16.11 shows that the shock S- (x) is a 5-entropy solution 
of (16.71) . which is obtained by finite-time blow-up as t — » CT from the analytic similarity 
solution (16.91) . 



2 In this connection, the result (ii) in Proposition ^. 2l sounds unusual: for non-analytic and very singular 
data, there exists a global analytic solution. 
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Shocks for non- degenerate NDE. For the corresponding non-degenerate NDE 
(6.18) u m = ((1 + u 2 )u x ) xx , 



the similarity solutions (16. 9p lead, on double integration, to the ODE (cf. (16. lip ) 

Az+B 



(6-19) 9' = TT^ (A^ 



It is easy to show using the phase-plane, that for z = — > 1 (this gives a necessary extra 
singular point (z ,go) of the flow, where g$ = z% — 1), (16.191) admits analytic solutions 
g(z) satisfying g(±oo) = C± > with C_ > C+, so as t — > CT, we obtain the shock. 

6.3. Stationary entropy shocks for other higher-order in time NDEs. We now 

very briefly check entropy properties of the shocks S±(x) for the following NDEs of arbi- 
trary order: 

(6.20) D 2 t m+2 u = D 2 x m (uu x ) (m > 1). 
For m — 0, this gives the following simple NDE: 

(6.21) u tt = uu x , 

for which both shocks S± are obviously weak solutions, so one needs to identify which 
ones are entropy. Note that, as (16. 7p . the PDEs (I6.20p for any m > 1 obey the Cauchy- 
Kovalevskaya theorem, so a unique local semigroup of analytic solutions does exist. 

S-entropy 5_ via analytic TWs. For a change, we present 5-deformations by TWs 



u{x,t) = f x (x - At) =}► A 2m+2 / (2m+2) = \ (f 2 Y 2m+1 \ 



(6 ' 22) A— f' = -Ul-n => h(y) = ^^. 



2 

We then observe that 

(6.23) /a(?/) —* S-(y) as A — >• uniformly in R, 



or 



so that the stationary shock wave S-(x) is G-admissible and is <5-entropy, where the 
necessary 5-deformation is given by the TW (16.221) with A = 5. 

A similar (but not explicit) construction of 5-entropy solutions with convergence (16.231) 
is performed for other normal NDEs such as 

(6.24) D 2m+i u = -D 2 x m {uu x ), or D 2m+2k u = (-l) k+1 D 2 x m {uu x ), k > 1. 

The corresponding analytic TW profiles f\(y) satisfying the convergence (16.231) in Ll oc 
are described in [T3J § 4]. 

Remark: S + can be formally created by a classical but non-analytic blow-up 
self-similar solution. There exists a self-similar blow-up to S + for the NDE (16.2ip via 



z > 

(6.25) u4x,t)=g(z), z = ^ =► g'z 2 = \ g 2 , so g(z) = \ ^ " ' 

[ l-4z' Z — U ' 
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This ODE obeys the symmetry 

(6.26) {» z :?; 

Note that g(z) is just C 1 (not C 2 ) at z = 0, which is enough to represent a weak solution of 
the degenerate PDE f)6.2ip (though, as we know, being weak often means almost nothing). 
Moreover, (16.251) is a classical solution of (I6.2ip . Observe that the non-analyticity of 
g(z) is associated with the too strong degeneracy at z = of the corresponding ordinary 
differential operator z 2 j-. We suspect that (I6.25j) is not entropy at all. Moreover, one 
can see that g(z) is not an odd function, so it looks more like a solution of an IB VP for 
x > with some boundary condition at x — 0. Nevertheless, we recall that it is a classical 
solution of the Cauchy problem. The NDEs (I6.20p deserve deeper study. 



7. On shocks for spatially higher-order NDEs 

7.1. Fifth-order NDEs. The similarity mechanism of shock formation remains valid for 
higher-order NDEs, among which, as an illustration, we comment on the following three 
(including the NDE-5 (051) ): 

(7.1) u a = —{uu x ) xxxx , 

l^ttt i^^x^xxxx- 

Concerning application of such fifth and higher-order NDEs, see [H], [TU p. 166], and 
references therein. The blow-up similarity solutions of S"_-type are the same, 

(7.2) u-(x,t)=g(z), z = x/{-tf, where f3 = ±, = §, = §, 
respectively. The ODEs are, respectively, 

(7.3) (gg ' ) W = -^ 7g ' z + 2 g''z 2 ), 

(gg')W = __!_(64 5 'Z + 57#V + V"2 3 ). 

These are much more complicated equations than all those studied before. We do not have 
a proof of existence of the 5*_-type profiles g(z) to say nothing about uniqueness, though 
we can justify that the shooting procedure to get a solution is well-posed according to 
dimensions of stable and unstable manifolds of orbits at the singular points z = (where 
g = 0) and z = — oo (where g = +1). On the other hand, the same numerical methods give 
us a strong evidence of existence-uniqueness. In Figure [T3| using bvp4c solver of MatLab, 
we present the unique solutions of the ODEs (17.31) satisfying the standard conditions 

(7.4) g{±oo) = =Fl and g(0) = g"{0) = # (4) (0) = (anti-symmetry). 

Note that first two ODEs admit solutions that are oscillatory about the equilibrium g = 1 
as z — > — oo, while the last one has monotone non-oscillatory solutions according to 
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Shock similarity profiles for u ,=..., u =..., u =-(uu ) and u ,=(uu ) 
3 r t tt ttt v x'xxxx t v x': 
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FIGURE 13. The shock similarity profiles g{z) for z < as the unique solutions 
of the problems ([7T5]> . ([731) and <E21)- 

the following asymptotics, respectively: for z <C —1, neglecting lower-order algebraic 
multipliers in the second and third formulae, 

g(z) - 1 ~ |z|-§cos (^5~^\z\i +c ), 

(7-5) g(z) - 1 ~ e-^ )§|z|f cos + c ), 

0(*)_l~e-§Cf)*W 4 . 

The exponentially small oscillations in the second line are hardly seen in the figure and 
requires another, logarithmic scale for revealing those. 

7.2. On a seventh-order NDE. For completeness and convenience of comparison, Fig- 
ure d3] also gives the shock similarity profiles (the dashed line) for the NDE-7, 

(7.6) ut = (uu x ) xxxxxx , where 

(7.7) u4x,t)=g(z), z = x/{-t) 1 7 => (gg')^ = ± g'z, g(±oo) = T l. 

The shock profile is very similar to that for the NDE-5 in (|1.15p . so that a general 
geometry of these shock profiles does not essentially depend on the order, (2m +1), of 
the PDEs ( I1.15P for m > 1; the oscillatory behaviour also changes slightly with m and 
always has the type given in the first line in (17. 5p . 

These results show that, for all the above higher-order NDEs, canonical shocks of S-- 
type are obtained by blow-up in finite time from smooth classical solutions. According to 
our 5-entropy approach, this confirms a correct entropy nature of such shock waves. 
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Shock similarity profiles for u=(uu ) ; a=-0.01 , 0, 0.01 , 0.02,. ..,0.05 
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FIGURE 14. Similarity profiles g(z) for z < satisfying the ODE in (|7.8p for 
various a £ [-0.01,0.05]. 

Let us describe other types of shocks and rarefaction waves for (17.61) driven by blow-up 
similarity patterns 

(7.8) u-(x,t) = (-t) a g(z), z = x/(-ty, (3=^ (gg') ^ = i±a </z - ag. 

These similarity profiles are presented in Figure [TU This shows that the profiles get more 
oscillatory for a < 0, but we failed to detect a "saw"-type profile as in [T8l § 4.3] for such 
a seventh-order ODE by using any numerical method. 

Finally, the analysis of the ODE in flTjjJ) on the invariant subspace (cf. the invariant 
subspace in [IBJ § 4.3]) W4 = Span{z, z 3 , z 5 , z 7 } shows that a nontrivial dynamics exists 
for the critical exponent 

^ = 11 = 0.161228..., 
and that the explicit solutions are given by 

g(z) = Cz + y|j z 7 , where C € R is arbitrary. 

8. On changing sign compactons for higher-order NDEs 

Finally, we return to the compacton solutions of the NDEs. First time, we discussed the 
entropy properties of compactons in Section H] for the NDE-3, where the entropy nature of 
such solutions was successfully justified. It turns out that the fact that these compactons 
are 5-entropy, i.e., are constructed by smooth ^-deformations, can be proved by a purely 
ODE approach, by smooth positive approximations of compactons via analytic solutions. 
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Compacton profile for u =-(u 2 u x ) xxxx +u 2 u x ; A.=-1/3, y. = 7 
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FIGURE 15. The first compacton profile of the ODE (pT2j) . 

We must admit that this ODE approach cannot be extended in principle to higher-order 
NDEs, so we need either to return to the original PDE 5-entropy method as in Section 
13. 2\ or to adapt the ODE approach to non-positive but less singular approximations (that 
we actually intend to do). 

8.1. Compacton for a cubic fifth-order NDE. For introducing a new model, un- 
like most of previous cases (excluding (11.111) in Section [575]) . without any hesitation, we 
consider the cubic NDE-5 

(8.1) u t = — {u 2 u x ) xxxx + u 2 u x in RxR + . 

We take the following TW compacton with the specially chosen wave speed A = — |: 

(8.2) u c (x,t) = f(y), y = x + \t -(/f > + f = / in R. 
We next perform the natural change leading to a simpler semilinear ODE, 

(8.3) F = f F (4) = F -F* in R. 

This easy looking equation admits a nontrivial countable set of various compactly sup- 
ported solutions that are analyzed by variational methods based on Lusternik-Schirel'man 
category and Pohozaev's fibering theory, [16]. Here we stress our attention to the primary 
facts that are connected with the proposed concepts of entropy solutions. 

The first and simplest compacton solution of the ODE (18. 2p is shown in Figure [TBI that 
was obtained numerically with the tolerances and regularization parameters 

Tols = lCT 10 and F^ ^ (z/ 2 + F 2 )~^F with also v = lCT 10 . 
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FIGURE 16. Enlarged zero structure from Figure [T5l near the right-hand interface. 

8.2. Oscillatory structure near interfaces: periodic orbits. In general, it looks 
that this compacton profile does not differ from those considered before as the explicit 
solutions in (14. 2 j) or (I5.10p . However, there is a fundamental difference that changes the 
mathematics of such solutions: for the fifth-order NDE (18. ip . the profiles f(y) = Fs(y) 
are oscillatory and are of changing sign near finite interfaces. In Figure [161 we show first 
three zeros near the interface at y = y$ > of the compacton profile from Figure [151 

In order to describe key features of such oscillatory behaviour at the right-hand inter- 
face, as y — > yo > 0, when F(y) — > 0, we perform an extra scaling by setting in the two 
leading terms of the ODE 

(8.4) F^=-F-s F{y) = {y -y) 6 <p(s), s = \n(y -y), 

where the oscillatory component ip(s) solves the following ODE: 

P A (ip) = e~ 2s [e- s (e- s (e- s (e- s (e 6 V) , ) , ) , ) , ] 
< ^ 8 ' 5 ' ) ee ipW + I6(f'" + 11V + 342^' + 36(V = -<^i 

It turns out that the oscillatory behaviour near the interface at y = y$ (i.e., at s = — oo) 
is given by a periodic solution v 9 *( s ) OI the ODE (18.51) . Namely, we list the following 
properties that lead to existence of a periodic orbit of changing sign: 

Proposition 8.1. The fourth- order dynamical system (\8.5§ satisfies: 

(i) no orbits are attracted to infinity as s +oo; 

(ii) it is a dissipative system with a bounded absorbing set; and 

(iii) a nontrivial periodic orbit (p*(s) exists. 

Proof, (i) The operator in (18.51) is asymptotically linear [2"2l p. 77] with the derivative at 
the point at infinity P4 that has the characteristic equation 

p 4 (A) = A 4 + 16A 3 + 119A 2 + 342A + 360 = (A + 6)(A + 5)(A + 4)(A + 3) = 0. 

35 



Therefore, all eigenvalues are real negative, —6, —5, —4, and —3, so infinity cannot attract 
orbits as s —>■ +oo. (ii) is a corollary of (i) after an extra scaling. 

(iii) Existence of a periodic orbit for dissipative systems is a standard result of degree 
theory; see [22], p. 235]. We complete the proof of existence by using a shooting argument 
as in § 7.1]. □ 

It turns out that the periodic solution y>*(s) is exponentially stable as s — > +oo (this is 
not easy to see from the ODE (18.51) by linearization and interpolation of the third term 
with the coefficient 119). The asymptotic stability of this periodic orbits is illustrated in 
figures in [T9J p. 187]. 

Thus, at the singular end point y = y^, the dynamical system (18. 3p generates a two- 
dimensional bundle of orbits with the behaviour 

(8.6) F(y) = {y- y Q f[ip*(s + s ) + ...], y e R+, s e R, 

where sq is an arbitrary phase shift of the periodic motion. Thus, the interface point 
{y = y$ , F = F' = F" = F'" = 0} is a complicated singular point (a zero) of the 
dynamical system f 1 8 . 3 j) . so one needs to check whether it corresponds to an entropy 
solution-compacton. It is worth mentioning that the 2D bundle ( 18. 6p matches with pre- 
cisely two symmetry conditions at the origin, F'(0) = F"'(0) = 0, and the existence of 
the compacton is confirmed by variational methods, [161 § 5]. 

8.3. Compactons are <5-entropy: a formal illustration. The oscillatory behaviour 
(18. 6p of the compacton near finite interfaces makes impossible to use the positive analytic 
^-approximation as for the NDE-3 in (14. 3p . Indeed, the same procedure for F = f 3 now 
leads to the "regularized" ODE 

(8.7) F s : F (4) = F - Fs + C s , where C s = 5^ - 5, so F s (y) -> 5 > 0, y -> oo. 

This gives the family {F$} consisting of functions Fg(y) that change sign finitely many 
times for all sufficiently small 5 > 0. These approximations Fs are less singular than the 
limit compacton profile F(y), which according to ( 18.61) is infinitely oscillatory as y — > y$ . 

The solvability of the approximating problem ( 18. 7ft can be traced out by the same 
variational method. Then the convergence 

(8.8) F s (y) -> F(y) as 5^0+ uniformly 

is associated with the stability of critical values of functionals; see [221 P- 387]. This 
^-approximation is shown in Figure UTT a). where the convergence (I8.8P is rather slow and 
is observed starting from 5 = 1CT 3 only, with the accuracy about 0.2. For 5 = 10~ 2 , the 
approximating profile F$(y) is still almost four times less than F(y) at the origin. The 
accuracy 0.1 is achieved for 6 = 10~ 5 . In (b), up to 5 = 10~ 8 , we show the zero structure 
of Fg(y) close to y , which, since F(y) ~ S > for y 1, is finite and each zero is 
transversal. These confirm that the approximating sequence {Fs}, though is of changing 
sign, is less singular than the compacton profile F(y) itself. 

An alternative approximating approach of such compactons is developed in [17], where 
F is approximated as e — ► + by the analytic family {F £ } of solutions of the regularized 
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ODE 

F £ : F (4) = F - (e 2 + F 2 )~^F (e > 0). 

Incidentally, this approach makes it possible to trace out the Sturmian index of some 
solutions by a homotopic connection to variational problems with known ordered set of 
critical points and known number of zeros for each of them, [IE] . 
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